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Abstract. Motivated by recent developments on solvable directed polymer 
models, we define a 'multi-layer' extension of the stochastic heat equation 
involving non-intersecting Brownian motions. 



1. Introduction 

We consider the stocliastic lieat equation in one dimension 

1 

with, initial condition Z{0,x,y) = S{x — y), where W{t,y) is space-time white 
noise (151 138) . The solution Z{t,x,y) is given by the chaos expansion 

Z{t,x,y) ^p{t,x,y) + 



(1) dtZ=-d^yZ + W{t,y)Z 



(2) y2 Pih,x,xi)p{t2-ti,xi,x2)---pit-tk,Xk,y) 

X W{dti,dxi) ■ ■ ■ W{dtk, dxk), 
where Ak{t) = {0 < h < ■ ■ ■ < tk < t} and 

For each t > and x,y G R, the expansion ^ is convergent in L2{W). It satisfies 
(IlJ in the sense that it satisfies the integral equation 

(3) Z{t,x,y)=p{t,x,y)+ I I p{t - s,y' ,y)Z{s,x,y')W{ds,dy'), 



which can be seen by iterating this equation. 

This solution arises as a scaling limit of discrete directed polymer models, 
in the 'intermediate disorder' regime [H I19j . It has also been known for some 
time that h = log Z arises as the scaling limit of the height profile of the weakly 
asymmetric simple exclusion process [6]. With this 'surface growth' interpretation. 
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h is understood to be the physically relevant solution (also known as the Cole-Hopf 
solution) to the KPZ equation |16j 



with 'narrow wedge' initial condition. 

In a remarkable recent development f2l[7l[l0l[lll[l2l[28l[29l[30l[3l]the exact 
distribution of Z{t, x, y) has been determined. The results so far have been based 
on two distinct approaches. One is to use the asymmetric simple exclusion process 
approximation together with recent work by Tracy and Widom j34L 1351 136L 137] in 
which exact formulas have been obtained for that process using the Bethe ansatz. 
The other is based on replicas, where the moments of the partition function are re- 
lated to the attractive 5-Bose gas and also computed using the Bethe ansatz. These 
developments indicate that there is an underlying 'integrable' structure behind KPZ 
and the stochastic heat equation which is not yet fully understood. 

It has recently been found that there exist exactly solvable discrete directed 
polymer models [HI HOI [22l [23l [32l [33] , yielding yet another approach. We will 
describe one of the main results from [22] . which provides the motivation for the 
present work. Define an 'up/right path' in R x Z to be an increasing path which 
either proceeds to the right or jumps up by one unit. For each sequence < si < 
• • • < SAr_i < t we can associate an up/right path (j) from (0, 1) to {t, N) which has 
jumps between the points (si, i) and {si, i+l), for i = 1, . . . , iV— 1, and is continuous 
otherwise. Let B{t) = {Bi{t), . . . , _Bjv(t)), t > 0, be a standard Brownian motion 
in and define 



and the integral is with respect to Lebesgue measure on the Euclidean set of all 
such paths. This is the partition function for the model. In |22| a formula is given 
for the Laplace transform of the distribution of Z^{t), which is obtained via the 
following 'multi-layer' construction. For n = 1, 2, . . . , iV, define 



where the integral is with respect to Lebesgue measure on the Euclidean set of 
n-tuples of non- intersecting (disjoint) up/right paths with respective initial points 
(0, 1), . . . , (0, n) and respective end points {t,N — n + l), . . . ,{t, N). Define {t) = 



logZf (t) and, for n > 2, {t) = \og[Z^ {t) / Z^_^{t)]. The relevance of this 



construction is analogous to the role of the RSK correspondence in the study of 
last passage percolation and longest increasing subsequence problems; in this setting 
it is based on a geometric (or 'tropical') variant of the RSK correspondence. The 
main resuh in [22 is that the process X^ {t) = {X^{t), . . . ,X^{t)), t > 0, is a 
diffusion process in M.^ with infinitesimal generator given by 



dth^-dlh+-idyhf + W{t,y), 




where 



E{(f>) = Bi(si) + B2{s2) - S2(si) + • • • + B^it) - Bn{sn-i) 




-A + Vlogi^o • V 
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where ipo{x) is the ground state eigenfunction of the quantum Toda lattice Hamil- 
tonian 

Af-l 

2=1 

The law of the partition function Z^(t) is obtained as a corollary. This result 
has been extended to a discrete-time framework in [9], where it is related to the 
solvability of a lattice directed polymer model with log-gamma weights introduced 
by Seppalainen [32j , also involves the eigenfunctions of the quantum Toda lattice 
(known as Whittaker functions) and works directly in the setting of the 'tropical 
RSK correspondence' introduced and studied in the papers |18L I21j . These results 
suggest that the continuum versions of the partition functions Z^{t), which we 
introduce in this paper, will play a role in our understanding of the 'integrable' 
structure which appears to lie behind KPZ and the stochastic heat equation. 

The outline of the paper is as follows. In the next section, we define the contin- 
uum versions of the above partition functions as chaos expansions involving non- 
intersecting Brownian motions. In section 3, we study the analogue of the partition 
functions when the space-time white noise is replaced by a smooth time-varying 
potential. In this setting we establish a connection with Darboux transformations 
of solutions to the heat equation, which give rise to evolution equations for the 
multi-layer process of partition functions. These equations are not directly mean- 
ingful in the white noise setting, but suggest that the multi-layer process has a 
Markovian evolution. There is a multi-dimensional version of the stochastic heat 
equation based on Brownian motion in a Weyl chamber, for which a Markovian 
evolution is readily apparent, and in section 4, we establish an analogue of the 
Karlin-McGregor formula for the solution of this equation. In section 5 we explain 
how, using this formula, the multi-dimensional process can (in principle — we only 
prove it for the first two layers) be expressed in terms of the multi-layer process, 
thus explaining the Markov property of the latter. We conclude in section 6 by 
formulating a continuum version of the tropical RSK correspondence. 

Acknowledgements. We thank Ivan Corwin, Martin Hairer, Jeremy Quastel, 
Gregorio Moreno and Roger Tribe for several helpful discussions. This research is 
supported in part by EPSRC grant EP/I014829/1. 



2. The continuum partition functions 



In this section we define the continuum analogues of the partition functions 
described in the introduction. For n — 1,2, . . ., t > and x, y eM., define 

(4) Zn{t,x,y)^p{t,x,y)"(l + f] f f R^^\{t^,x^), . . . ,{tk,Xk)) 

X W{dti,dxi)-- -Width, dxk)^, 

in) 

where i?^, is the fc-point correlation function for a collection of n non-intersecting 
Brownian bridges which all start at x at time and all end at y at time t. Note 
that 2,1 = Z is the solution of the stochastic heat equation defined by (|2|) . 

Theorem 2.1. The series ^ is convergent in L2{W). 



4 



NEIL O'CONNELL AND JON WARREN 



Proof. Wc need to show that 



(5) 



y-f I 



R^^'{{ti,xi), {tk,Xk))'^dxi ■ ■ ■ dxkdti ■■ -dtk <oo. 



This is equivalent to showing that Ee^ < oo, where L is the total intersection local 
time between two independent copies X = < s < t, i = and 

Y = {Yg, < s < i = 1, . . . ,n) of the system of n non-intersecting Brownian 
bridges which all start at x at time and all end at y at time t. We will show 
that, in fact, all exponential moments of L are finite. Without loss of generality 
we can assume that x = y. First note that L = A + B, where A is the intersection 
local time on the time interval [0,t/2] and B is the remainder; by symmetry, A 
and B have the same distribution. Thus, by Cauchy-Schwartz, it suffices to show 
that A has finite exponential moments of all orders. Now, on the time interval 
[0, t/2], X and Y are (up to a bounded Radon-Nikodym density and time-change) 
two independent copies of Dyson Brownian motion in A„ started at i. It therefore 
suffices to show that, for two independent Dyson Brownian motions in A started 
at the origin and run for time T, say, the total intersection local time has finite 
exponential moments of all ordcirs. Dc^notc; these two Dyson Brownian motions by 
U and V. Although these processes are necessarily defined via an entrance law (as 
they are started on the boundary of the Weyl chamber), it can be shown that they 
satisfy a system of SDEs 

du: = di3i + j2-^, dv: = drd+y- 



v: - vi 



where /3*, 7% i = 1, 2, . . . , n are a collection of independent standard one-dimensional 
Brownian motions. The total intersection local time up to time T is given by 
Y^i^jL^^, where denotes the local time (at zero) of — up to time T. 
Again by Cauchy-Schwartz, it suffices to show that for each distinct pair 
has finite exponential moments of all orders. In the following we will use the fact 
that the random variables \U^\ (for each i) and the absolute value of a Gaussian 
random variable all have finite exponential moments of all orders. 

By Tanaka's formula, 

= |C/^ - V^\ - r sgn(C/j - Vi)d{Ui - Vi) 

Jo 

= m-Vi\- [ sgn(/7:-y/)rf(/3:-7i)- / sgn{U:-V,^){Di-Ei)ds 
Jo Jo 



sgn([/:-y/)d(/3:-7^) 



^\Dl\ds+ r \Ei\ds, 
Jo 



where 



Vi - Vi 



Thus, it sufiices to show that each of the random variables 



m-vi\, 



f sgn{Ul- 
Jo 



VJ)d{Pl - il) 
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/ \Di\ds and / |E^|ds, 
Jo Jo 
has finite exponential moments of all orders. For the first, we note that 

\ui.-v^\ < \ui,\ + \v^\, 

which has finite exponential moments of all orders. The second is the absolute value 
of a Gaussian random variable with mean zero and variance 2T. So it remains to 
show that, for each i, 

rT 



■■= I \Di\ds 
Jo 



has finite exponential moments of all orders. We will prove this by induction over 
i. For i < j, define 

= lo wfhji''- 

First we note that 



Ci = / Dlds = U^- p}r, 
Jo 

has finite exponential moments of all orders. Now, since 

a = Ci2 + • • • + Cm 



and each term is non-negative, this implies that ^ij < ^ and hence has finite 
exponential moments of all orders for each j = 2, . . . ,n. Now 

' — + ^i;l2• 



6 = 62 + 63 + • • • + ^2n = I Dlds + 2^12 = U^-I5l + 2ii 







Thus ^2 and ^23, • • • , $2n all have finite exponential moments of all orders. Similarly, 

6 = 63 + 63 + 64 + • • • + 6n = / Dlds + 2^3 + 2C23 = C/| - + 2^13 + 2^3, 



and so on. □ 



3. Darboux transformations and non-intersecting Browniain motions 

In this section we replace the white noise potential by a smooth potential (j), 
which wc assume for convenience to be in the Schwartz space E of rapidly decreasing 
smooth {C°°) functions on M+ x K. 

For each n = l,2, ...,t>0 and x, j/ G R, define 

(6) Zn{t,x,y)=p{t,x,yr¥.e^v (^J^ cl>{s,Xi)ds^ , 

where XI, < s < t, i = 1, . . . ,n denote the trajectories of n non-intersecting 
Brownian bridges which all start at x at time and all end at y at time t. On 
one hand, these are the analogues of the partition functions Z„'s introdiiccd in the 
previous section with the white noise replaced by a smooth potential. On the other, 
they are directly related to the Z„'s by the formula 

(7) Zn{t, x,y)=E [Zn{t, X, y) exp<^(W^(<^))] , 
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where exp*(VK(0)) is the Wick exponential of W{(f>) defined by 

exp^(M^((/))) = exp ( /" / (l){s,x)W(ds,dx) j [ (p{s,xfdxds 

In other words, as a function of 4>, Zn{t, x, y) is the S'-transform of x, y) 
To see that ^ holds, on the RHS replace x, y) by the series (|4]) and exp*(W^((/))) 
by its Wiener chaos expansion; computing the expectation of the product of these 
two series we obtain 

p{t,x,yY^^ / / (j}{ti,xi)...(l){tk,Xk)R^^\{ti,Xi),...{tk,Xn))dxi...dxkdti.. 
i^^QJA^t) Jr." 

= p{t,x,yrEc^j> (j^l <l^i^,Xl)ds\ . 



By the Feyman-Kac formula, Z := Z\ satisfies the heat equation 

(8) dtZ=h)lZ^^{t,y)Z 
with initial condition ^(0, a;, = b(x — y). 
Proposition 3.1. 

(9) Z„(t, X, y) = cn,t det [didlZ{t, x, y)] , 
where c„,t = i"(»-i)/2 l^ll j[. 



We will prove this via a generalisation of the Karlin-McGregor formula. Set 

(10) A„ = {x e R" : xi > • • • > x„}. 
For each t > and x, y € A° , define 

(11) Z„(i,x,2/) =p*(t,x,y)Eexp ^^y" (/)(s, C/*)ds^ , 

where ?7 is a collection of non-intersecting Brownian bridges started at positions 
xi, . . . , x„ and ending at yi, . . . , yn, and Pn(t, x, y) is the transition density of a 
Brownian motion in A„ killed when it hits the boundary, given by the Karlin- 
McGregor formula jl7j . 

n 

P*n(t,x,y)^ ^ sgn(cr) J|p(t,x,,y^(,)). 

Proposition 3.2. 

(12) Z„(<,x,y) =det[Z(i,x„2/,)]^j-^i. 

Proof. According to the Feynman-Kac formula, Z„ satisfies the equation 

(13) dtZn = \/^yZn + J2 '^(^' y^)^" 

with Dirichlet boundary conditions on 9A„ and initial condition Z„(0,x,y) = 
(S(xi — yi). Moreover it is the unique solution to this initial-boundary value 
problem which vanishes as |y| — >■ oo uniformly for t in compact intervals. 
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On the other hand, dct[Z{t, Xi,yj)]fj^i satisfies the same initial-boundary value 
problem and vanishes as \y\ — >■ co uniformly for t in compact intervals. So the 
identity follows by uniqueness. □ 



Proof of Proposition 13.11 For a e M, denote by a e M" the vector with 
all coordinates equal to a. Now, it is immediate from the definitions that 

Zn{t,a,b) Zn{t,x,y) 
— -— — lim -. 

By Proposition 13 . 21 we have 



Zn{t,x,y) 
hm . , ... . — det 
x^a,y^b A[x)A(y) 



dldiZ{t,a,b) 



n-l 



i,j=0 

where A{x) = ni<j(^i ~ ^j) (^^6, for example, [39]). On the other hand, 
x^a.y^b A{x)A(y) 

which completes the proof. □ 

Define Un{t, x, y) recursively by Z„ — uiU2 ■ ■ • Un- 
Proposition 3.3. The functions u„ satisfy the coupled system of heat equations 

(14) dtun = \dlu^ + m, y) + dl log (Z„_i/p"-i)]u„ 

with initial conditions Un{0,x,y) = S{x — y). 



Proof. The equations follow from Proposition [32] together with known prop- 
erties of Darboux transformations of solutions to one-dimensional heat equations 
with time- varying potentials, see for example [3]. The initial condition follows 
immediately from the definition of Z„. □ 



The coupled heat equations of Proposition 13731 are not immediately meaningful 
if we replace the smooth potential 4> by space-time white noise. However they do 
suggest that the multi-layer process 

a;, •),••■, Z„(i,x,-)), t>0 

is Markov. In the following, we introduce a natural extension of the multi-layer 
process which will play an important role in our understanding of the Markov 
property when we return to the white noise setting. 

Define, for t > and x, ?/ G A„, 

(15) ^"(*'^'2')= A(x)A(y)- 

This extends continuously to the boundary of A„ x A„; by Proposition 13.21 for 

x€R, 

(16) Znit,x,y) = A(2/)-Met [dl-'Z{t,x,y,)]l^^^ . 

Rather surprisingly, we will now show that the apparently richer object Zn{t,x, •) 
is, for a fixed x G R and t > 0, given as a function of {Zi(t, x, •),... , Zn{t, x, •)). 
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For notational convenience, set Zq = 1 and Ai = M. For n > 1, define 

nt Zl 

We note that Sn — dxy log Zn, from Lemma IXHl below. For z £ A„_i and y e A„, 
write z ^ y ii yi > zi > 1/2 > ■ ■ ■ > yn~i > Zn-i > yn- For y G A° , denote by 
GT{y) the Gelfand-Tsetlin polytope 

{{y\y\ . . . , y"-^) e Ai x A2 x • • • x A„_i : y^ ^ y^ ^ . . . ^ y'-i -< y}. 

Theorem 3.4. For t > 0, x eR and y e A° , 

n n— 1 n — 



In ttie case ip — 0, tliis reduces to tlie fact that the volume of GT{y) is pro- 
portional to A(?/). By pB)) and Proposition 13.11 this theorem can be seen as a 
consequence of the next two lemmas and the fact that Cn-i,tCn+i,t/ cl^ t — '^^i 

Lemma 3.5. ///i, /2, ■ • • is a sequence of continuously dijferentiable functions on 
R with fi = l then 

demy,)]ls-i = / det[/;+i(z,)]i:7iid^i ■ • • rf^n-i- 

Jz^y 



Proof. Using the formula 



Iz-iy - dot [lyj + i<2:i<iyj j j^i 



we have, by Cauchy-Binct, 

det[/,'+i(zj)]"7iidzi . . . 



z^y 



1 n-1 



det[/;+i(zj)]-j-ii det [h,+i<z^<yAljli '^^i ■ ■ • '^^"-1 



f[+i{z)dz 



det 

det[/.(2;,)]^,^i. 



det [f;+i{yj) - /»+i(2/j+i)]-j-ii 



as required. 



□ 



Lemma 3.6. Let g(x,y) be a smooth function and define Wq — I, Wi = g and, 
for n>2, Wn = det [9* 9^(7(0;, y)] "^.j'^^. Suppose that Wn is strictly positive for all 
n > 1 and define Tn = Wn-iWn+i/Wf^. Then the following identities hold: 

T„ - d^ylOgWn = dy{dy{- ■ ■ dyid^ g / / T,) / T^) ■ ■ ■ )/T„_l), 



71 — 1 71 — k 



I L n I L ± I L ri 

det[dl''gix, y,)]l^^, = [] g{x, y,) [] 11 ^^(^' V^'"')'^ 



1 1=1 
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Proof. From well-known properties of Wronskian determinants, dxWn, dyWn 
and dxyWn can be expressed as determinants, namely 

dxWn = det j=o,...,n-i ' 

dyWn = det [5;5^5(a;,2/)],^o,...,„_i^ ,=o,i.....n-2,„ ' 
dxv^n = det\dldlg(x,y)] . ^ . ^ . 

f L ^ ^ ' "'Ji=0,l,...,n-2,n; j=0,l,...,n — 2,n 

It follows from Sylvester's determinant identity |141 p22] that 

WndxyWn - idxWn)idyWn) = Wn^lWn+1, 

proving the first identity. Similarly, for any fc > 1, 



This implies that 



In particular. 



{dy{d:g/g))/T, = {d:-'W2)/W2, 
{dy{d:-'W2/W2))lT^ = (9r'M^3)/W^3 



and so on, yielding the second identity. 
Now, by Lemma 13.51 



det 



det 



n-l 



n ^r' 



=1 i=l 



det 



y" ^<y 



yj g{x,y"~^) 
Ti(x,yr') 



n~l 



i=l 



Applying Lemma 15751 again, using Ti — dy{dxg/g), we obtain 



det 



(9,-1— ^^TT— 

g{x:yj ) 
Ti{x,yr') 



det 



d. 



g{x,yj ) 



n-2 



n-2 



det 



aj+lg(x,B""^) 

o 

T2ix,yr') 



n-2 



n-2 



l[T,{x,yr')dyr 



Now apply Lemma 1331 again, using T2 = dy{dy{d^g/ g)/Ti), and so on, to obtain 
the third identity. □ 

The evolution of the Sn's is given by the following proposition. 
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Proposition 3.7. For n > 1, 

(17) dtSn = Id^Sn + dy [Sndy log Ur.] . 



Proof. From Proposition 13.31 the functions hn ~ logu,i satisfy 

dthn = \dlK + \{dyKf + cf>{t, y) + log (Z„_i/p"-i) . 
It follows immediately that 5*1 = dxyhi satisfies 

dtSl = ^d^Sl+dy[Sldyhl]. 

We prove the general statement by induction. Assume the induction hypothesis 

dtSn = ^d^Sn+dy[Sndyhn]. 



Now Sn+i = Sn + dxyhn+i and 

-dyhn+i 4 -y^y 



dtK+i = IdlK+i + \{dyK+if + ^{t, y) + dl log (Z„/p") . 



Thus, 

dtSn+l = T,dySn+l + dy[Sndyhn] + dy[{Sn+l ~ Sn)dyhn+\\ + dySn- 



But 

dySn = dy[Sndy \0gSn] = dy[Sn[dyhn+i - Oyhn)], 

SO we have ^ 

dtSn+l — —dySn+i + dy[Sn+ldyhn^i], 

as required. □ 



4. The Karlin-McGregor formula 



For n — 1,2, . . . and x,y £ A° , define 
(18) Zn{t,x,y)^p*Jt,x,y)(l + Y, / ((^i, xi), . . . , (tfe, x^)) 

X Widh,dxi)---W{dtk,dxk)y 

where i?^"'' is the fc-point correlation function for a collection of n non-intersecting 
Brownian bridges started at positions x and ending at positions y at time t. 

Proposition 4.1. The series il8]} is convergent in L^{W). 

Proof. We need to show that E[e'^l/i] < oo where L is the total intersection 
local time between two independent sets of n independent Brownian bridges started 
at positions x and ending at positions y at time t, and A is the event that each 
set is non- intersecting. So it suffices to show that Ee'^ < oo. By considering 
pairwise intersection local times and applying Holder's inequality one obtains Ee^ < 

Ee"''V*72^ where R is the local time at zero of a standard Brownian bridge on 
[0,1], which has the Rayleigh distribution P(i? > r) = e""" r > (see, for 
example, [24]). □ 
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For each n, the function Zn{t,x,y) satisfies the equation 

i 

with initial condition Z„(0, x, y) — J^. 5{xi — yi) and Dirichlet boundary conditions 

on 9A„. By this we mean that Z.„(t, x, y) satisfies the integral equation 

(19) 

Z„(t,a;,y) V / / pl{t ~ s,y' ,y)Z^{s,x,y')W{ds,dy'^)\\dyr. 

Iterating the integral equation (|19l) and then integrating out the y'^ variables which 
don't appear in the white noise yields the chaos expansion Moreover, for each 
X € A°, X, •), i > is a Markov process, which can be seen as a consequence 
of the 'flow property' 



Zn{s + t,x,y) = / Zn{s,x,z)Znit,z,y;s), 

"'A„ 

where 2.„(t, x, y; s) is defined via the chaos expansion (jlSp but with the shifted 
white noise W{s + •, •). This flow property can be seen directly from the chaos 
expansion or by taking S'-transforms of both sides and using independence. 

Theorem 4.2. For x,y e A°, Zn(t,x,y) = dct[Z(t,Xi,yj)]lj^j^. 

Proof. Let (f) £ E, multiply both sides by exp^ W{(j)) and take expectations. 
The LHS becomes Zn{t, x, y), deflned earlier by pTt . which satisfies 



dtZn = ^\Zn + ^ (t){t, yi)Zn, 



with initial condition Zn(0, x, y) = 5[xi — yi) and Dirichlet boundary conditions 
Zn{t, a;, y) = for i > and y e 9A„. The RHS becomes 

Cnit, X, y) := E [det[Z(<, a;., yj)]"j=i exp*(iy(0))] , 

which is given as follows. Let . . . be a collection of independent Brownian 
bridges which start at positions x and end at positions y at time t. Define 

Vn{t,x,y) =p„(t,x,2/)Ecxp j (jjis^Pl) + L* j , 

where p„{t,x,y) = JJiPit, Xi,yi) and 



denotes the total intersection local time between the Brownian bridges. Note that 
Vn{t,x,y) is defined for any x,?/ S R", and satisfies 

Vr,{t,x,y) lf[Z{t,x^,y,)exp^{W{^))] . 
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This can be seen by multiplying together the chaos expansions of each term in the 
product and taking expectations, which can be justified by the fact 5^ that each 
term is in LP{W) for all p > 1. It follows that 

(20) C^{t,x,y) = sgn{<j)Vn{t,<JX,y). 

Note that C„(0, x, y) = J| - 6{xi — yi) and, for i > 0, C„(t, x,y) —Q for y <E dKn. 
Now, by Feyman-Kac, for each cr G S'n, Vn{t,(7x,y) satisfies the equation 
dtu = ^AyU + ^ (j){t, yi)u + ^ 5{yi - yj)u, 

i i<j 

in its integral form; more precisely, 

Vn{t,ax,y) = pn{t,ax,y) + Y / / Pnit - s,y' ,y)Vn{s,ax,y')(j){s,y'i)dy'ds 

^ Jo JR" 

+ X] / / Pn{i- s,y' ,y)Vn{s,(Jx,y')5{y[-y'j)dy'ds. 

i<j-^0 ^R" 

By linearity, C„ (i, x, y) (defined by (PO)) for y e M") also solves this integral equation 
and by the Dirichlet boundary conditions satisfied by C„(t, x, y) we obtain 

Cn{t,x,y)^Pn{t,x,y) + y2 I / Pn{t- s,y' ,y)Cn{s,x,y')(l){s,y[)dy'ds 

■ Jo JR" 

+ / / Pnit- s,y',y)Cn{s,x,y')6{yi~y'^)dy'ds 

»<j-^0 ^R" 

= Pnit,x,y) + Y / Pn{t- s,y' ,y)Cn{s,x,y')(l){s,y[)dy'ds. 

■ Jo jR" 

It follows that Cn{t,x,y) satisfies 

with initial condition C„(0, a;, y) = J^^ S{xi — yi) and Dirichlet boundary conditions 
on 9A„. 

Now we will argue that Vn{t,x,y), and hence Cn{t,x,y), vanishes as \y\ — >■ oo 
uniformly for t in any compact interval. Indeed, since is bounded, we have for 

t e [0,T] say, 

|K(i,x,2/)| <p„(i,x,y)e^"^Ee^*, 
where C is a constant. It is straightforward to obtain a bound on Eexp(Lt) which 
is uniform in t € [0,T] and x,y G R": by considering pairwise intersection local 
times and applying Holder's inequality one obtains 

where R is the local time at zero of a standard Brownian bridge on [0, 1]. Thus 
Vn{t,x,y) — >■ as |?/| — oo uniformly for t e [0,T], as required. 

By uniqueness, we conclude that Cn{t,x,y) — Zn{t,x,y) for < > and x,y £ 
A° . Since this holds for any (j) G E, we are done. □ 



A MULTI-LAYER EXTENSION OF THE STOCHASTIC HEAT EQUATION 



13 



5. On the evolution of the partition functions 



In this section we discuss the analogue of Theoreni l3.4l in the white noise setting, 
and the imphcation that {Zi(t, a;, •),... , Z„(t, x, •)), i > is Markov. 

We expect, but will not prove here, that for each i > 0, 

has a version which almost surely extends continuously to a strictly positive function 
on A„ X A„. In particular, for each t > 0, almost surely, 

(21) Zn{t,a,b) ^ Cn,t lim Z„(t, a:, y), 

uniformly on compact intervals. Assuming this continuity it can be shown that the 
analogue of Theorem 13.41 holds in the white-noise setting, that is, if we set Zq = 1 
and define, for n > 1, 

p f '^n— l^n+l 

" nt Zl ' 
then, for t > 0, x G M and y G A° , 

n „ n—ln—k 

(22) Z^{t,x,y) = ^{y)-^\{ZS,x,y{) X{X{^kit,x,yr'')dyr'' ■ 



fe=l i=l 



It is not difficult to sec that for each a; G M and for each n, the process Z„(t, x, •), 
t > has the Markov property. Assuming the validity the formulas ([21]) and (|22|) 
this would imply that {Zi{t, x, •),... , Z„(t, x, •)), t > is a Markov process. 

A proof of the existence of an almost surely continuous extension for Z„(t, x, y) 
based on Kolmogorov's criterion would be long and technical. Here we will satisfy 
ourselves with a continuous extension in L2, which then allows us to prove ([22l) . 
and hence the Markov property of the multi-layer process, in the special case n = 2. 

Lemma 5.1. For each t > 0, 

A / , N _ '^njt, X, y) 

A{x)A{y) 

extends continuously in L2{W) to A„ x A„. Moreover this extension satisfies 

Zn{t,x,y) = Cn,tZn{t,X,y). 

Proof. First we recall that we have the representation 
(23) Ut,x,y) = f.^Y'^\ (l + f2 f f Rt\{ti,x,\...Atk,xu)) 



X W{dti,dxi) ■ ■ ■ W{dtk,dxk) 



~ (n) 

where Rj. are the correlations functions of a collection of n non-inersecting Brow- 
nian bridges starting at x and ending at time t at y. Since 

Pn{t,x,y) 
A{x)A{y) 
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extends continuously to A„ x A„ this representation naturally defines the extension 
of 2, to A„ X A„. Our task is show continuity in L2{W). For this it is enough to 
show that 

{x,y,x',y') E[Zn{t,x,y)Zr,{t,x' ,y')] 
is continuous. Now, similarly to as in Theorem 13.41 this expectation is equal to 

' pUt,x,y) y ^ 
A(-)A(y)j "^^^ ^ 

where L is the total intersection local time of two independent sets of non-intersection 
bridges, X and X' say. 

Let us write L = L[o + L^s i_g-^ + L^t.^ ^j, where L[o,e] denotes the local time 
accrued over the times periods [0,(5] and so on. By conditioning on the position of 
the bridges at times S and t — S we have 

E[cMLis,t-5])\iX{0),X'{0),X{t),X'{t)) = z] = 

p{z, C)E[exp(L[5,t„,])|(X(,5), X'{5), X{t - 6),X'{t - 5)) = C] dQ. 

The kernel p{z, C) can be written as a product of transition densities for non- 
intersecting Brownian motions, and is thus seen to be continuous. From this it 
follows by a dominated convergence argument that 

E[eML[s.t-s])\{X(Q),X'{0),X(t),X'(i)) = z] 

depends continuously on z also. 

To deduce the continuity of z E[exp(i:)|(A:(0), X'(0), A"(t), ^ z] we 

must show that the difference 

E[exp(L)|(X(0),X'(0),X(t),X'(i)) = z] 

- E[exp(L[5,t_5])l (^(0), X'(0), X'{t)) = z] 

is can be made uniformly small for z within compact sets by choosing 5 small 
enough. Applying the Cauchy-Schwartz inequality this amounts to showing that 

E[exp(4i[o,5])l(^(0),^'(0),X(i),X'(t)) = z] 

and 

E[exp(L[t_5^,])|(X(0),X'(0),X(t),X'(i)) = z] 

can be made uniformly close to 1, and this can be achieved using similar calculations 
to those made in Theoreml3.4l □ 



Theorem 5.2. For x e M and y e A2, 

^24) Z2{t,x,y) ^ 1 Z2{t,x,z) ^_^ 

Z{x,yi)Z{x,y2) yi-y2jy^ Z{t,x,zy 

Proof. It is known |5j that the solution to the stochastic equation Z{t,x,y) 
admits a version that is almost surely continuous in t and y and moreover strictly 
positive. We assume in the following that we are using this version. In particular, 
having fixed t, x and yi > y2 we let Ae{x) be the event {Z{t,x,z) > e for all z e 
[2/2, yi + 1]}. Then as e J, we have P(A,(a;)) f 1- 
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By Theorem 1121 for u,v e A^, 
1 



Z2(t, U, V 



Hence 

(25) 



' A{u)A{v) 

Z2{t,U,v) 



[Z{t,Ui,Vi)Z{t,U2,V2) - Z{t,Ui,V2)Z{t,U2,Vi)]. 

Z{t,Ui,Vi) Z(t, Ui,W2)" 



Z{t,U2,Vl) Z(t,U2,V2) 



Z(t,u2,vi)Z{t,u2,v2) A{u)A{v) 
Let h > 0, V = {z + h, z) and integrate this equation with respect to z over the 
interval [t/2,2/i] thus obtaining the identity 



Z2{t,u,{z + h,z)) 

Z(t, U2, Z + h)Z(t, U2, z) 



dz 



1 



(ui - U2)h 



y^+''z{t,u,,z)^ fy-+'' z{t,u^,z) 



Z{t,U2, Z 



■dz 



Z(t,U2,z) 



dz 



Now let h tend to zero. By the continuity of Z{ui, •) and Z{u2, •) the RHS converges 
almost surely to 



1 



(ui - U2) 



Zit,ui,yi) Z{t,ui,y2) 



Z{t,U2,yi) Z(t,u2,y2) 



We want to identify the limit of the LHS. Consider 



E 



-dz 



< 



Z2{t,U,{z + h,z)) 
Z{t,U2,Z + h)Z{t,U2,z) 

y^ \Z2{t, u, {z + h, z)) - Z2{t, u, z)\ 



y^ Z2{t,u,z) 

Z{t,U2,ZY 

dz 



dz 



y2 Z{t,U2,Zy 
y^ \Z{t, U2,Z + h) — Z{t, U2, z)\Z2{t, u, (z + h, z)) 



S/2 



Z{t, U2, zyZ{t, U2, z + h) 



dz 



We have 



(26) E[E; A,{u2)] < e^^ E\Z2{t,u,{z + h,z)) - Z2{t,u,z)\dz 

{E[Z2{t, u, {z + h, z)f]E[{Z{t, U2, z + h)- Z{t, ua, z)f])^^'^)dz 



By virtue of the uniform continuity in L2 of the mappings (zi, Z2) '-^ Z2{t, u, (21, Z2)) 
and z I— > Z(t, ui, z) these integrals tends to zero as, h and consequently E tends 
to in probability. Thus we have proven 



(27) 



1 



(ui - U2) 



Z{t,ui,yi) Z(t, ui,?/2) 



Z2it,u,z)) 
Z{t, U2, zY 



dz. 



Z{t,u2,y\) Z{t,u2,y2)_ 

Next let u = {x + h,x) and let h 10. The LHS of ((27)) can be rewritten as 

_ Z{t,{x + h,x),{yi,y2)) _ 
Z{t,x + h,yi)Z{t,x,y2) ' 
as h \.Q this converges in probability to 

Z{t,x, (2/1, ya)) 



(yi - y2] 



Z{t,x,yi)Z{t,x,y2)' 
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On the other hand, if we consider 
F = 

we have 



dz — I -— tttuz 



E[F;A(a;)] < e"^ / ^Z2{t,{x + h,x), z) - Z2{t,x, z)\dz 

which again by the L2 continuity of Z2 converges to as /i | 0. From this it fohows 
the RHS of (I22I) converges to 



z{t,x,zy 

and the result is proven. □ 



We remark that the identity p7p shows that the ratio of two solutions to the 
stochastic heat equation is in H^; in fact, it has recently been shown by Hairer |13| 
to be in C^/'^''. 

Corollary 5.3. For x E R, {Zi{t,x,-),Z2{t,x,-)), t > is a Markov process. 

It will also be interesting to understand the evolution of the multi-layer process 
in terms of a system of SPDEs. Motivated by the evolution equations obtained in 
Section 3 in the the case of a smooth potential, it is natural to consider and to try 
to make sense of the system of equations 

(28) dtBn = \d^Sn + dy [§ndy log U„] , 

where Zn — Hi ■ ■ ■ ILn ■ For recent progress in this direction, see |13| . 

6. Concluding remarks 

As remarked in the introduction, the multi-layer construction presented in this 
paper is based on the tropical RSK correspondence, so it is natural to consider such 
an interpretation in the continuum setting. The analogue of RSK in the context of 
section 3 is the mapping ^^^^ h-> {u„(t, 0, •), n > 1}. In the language of RSK, 
{un{t, 0,x), n > 1; a; > 0} is the P-tableau, {M„(t, 0, —x), n> 1; x > 0} is the Q- 
tableau, and their common 'shape' is the sequence {un{t, 0, 0), n > 1}. We note the 
following symmetry, which corresponds to a well known symmetry property of the 
RSK correspondence. Writing / = ^IjqjJxK' ^(/) ~ {wn(t, 0,x), n > 1; a; > 0}, 
Qif) = {un{t,0,—x), n > 1; X > 0} and /t(s, x) = f{s,—x), we have: P{f^) = 
Qif) and g(/t) P{f). 

Similarly, in the white noise setting, we define 

P(W^) = {Unit, 0, x), n > 1; x> 0} Q{W^) = {U„(t, 0, -x), n > 1; x > 0} 

where denotes the restriction of W to [0, t] x M. Given the main result of the 
paper [22j . it is natural to expect that P{W^) and Q{W*) are diffusion processes 
in (indexed by x > 0) which are conditionally independent given their starting 
position {![„(<, 0,0), n > 1}. This would be the analogue, in this setting, of 
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Pitman's '2M — theorem. In recent work, Corwin and Hammond (8j have 
shown that the process {'U„(i, 0, a;), n > 1; a; G K} has a kind of Gibbs property 
and, moreover, has all components strictly positive almost surely; they call the 
logarithm of this process the 'KPZ line ensemble'. For large t, it should rescale 
to the multi-layer Airy process (or 'Airy line ensemble') which was introduced by 
Prahofer and Spohn [25]. At present, we only know this to be the case for the one- 
point distributions of the first layer: it has been shown in the papers [2|, I31j that 
the distribution of \og[Z{t,0,x)/p{t,0,x)] (which is independent of x) converges 
in a suitable scaling limit to the Tracy- Widom distribution. This result on the 
first layer has been tentatively extended to the finite-dimenisonal distributions by 
Prolhac and Spohn [26l [27] . 
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